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Abstract 

The long awaited discovery of a new light scalar boson at the LHC opens up a new 
era of studies of the Higgs sector in the Standard Model and in its extensions. In this 
paper we discuss the consequences of the observation of a light Higgs boson with the mass 
and rates reported by the ATLAS and CMS collaborations on the parameter space of the 
phenomenological MSSM, including also the LHC searches for heavier Higgs bosons and 
supersymmetric particle partners, as well as the constraints from S-physics and dark 
matter. We explore the various regimes of the MSSM Higgs sector depending on the 
parameters Ma and tan /3 and show that only two of them are still allowed by all present 
experimental constraints: the decoupling regime in which there is only one light and 
standard-like Higgs boson while the heavier Higgs states decouple from gauge bosons, 
and the supersymmetric regime, in which there are light supersymmetric particle partners 
which might affect the decay properties of the light Higgs boson, in particular its di- 
photon and invisible decays. 



1 Introduction 



Results of the search for the Higgs bosons at the LHC with the 10 fb 1 data collected in 
2011 at 7 TeV and 2012 at 8 TeV have just been presented by the ATLAS [l] and CMS |2j 
collaborations and there is now a 5a evidence by each of the experiments for a new particle with 
a mass of ~ 126 GeV. Complementary evidence is also provided by the updated combination 
of the Higgs searches performed by the CDF and DO collaborations at the Tevatron [3] , which 
has also been just released. As we are entering an era of Higgs studies, these results have deep 
implications for the Minimal Supersymmetric Standard Model (MSSM). The implications of a 
Higgs boson with a mass value around 126 GeV in the context of Supersymmetry have been 
already widely discussed [4jj7] since the first evidence of a signal at the LHC was presented at 
the end of 2011. In particular, we have discussed the consequences of the value of M h for the 
unconstrained phenomenological MSSM (pMSSM) with 22 free parameters [8], for constrained 
MSSM scenarios such as the minimal gravity, gauge and anomaly mediated SUSY-breaking 
models, and in scenarios in which the supersymmetric spectrum is extremely heavy [2] . We have 
shown that only when the SUSY-breaking scale is very large or the mixing in the stop sector 
is significant the observed value can be accommodated in these models. This disfavours 
many constrained scenarios such as the minimal anomaly and gauge mediated SUSY-breaking 
models and some (even more constrained) versions of the minimal super-gravity model. 

In this paper, we extend the previous study in new directions afforded by the improved 
data from the LHC experiments. First, we refine our analysis of the implications of the value 
Mh ^ 126 GeV for the decoupling regime by considering different types of stop mixing scenarios 
which significantly affect the maximal mass value and we explore the implications of a broader 
range of the top quark mass value, rrit = (173±3) GeV, on Mh- Then, we analyse in detail the 
implications of the ATLAS and CMS searches for the heavier MSSM Higgs bosons, the CP-even 
H, a pseudoscalar A and two charged H ± states. In particular, we discuss the A/H/h — > t + t~ 
for the neutral and the t — > bH + — > bru searches for the charged states to further constrain 
the [Ma, tan /3] parameter space, where tan/3 is the ratio of the vacuum expectation values 
of the two Higgs doublet fields. We also discuss the effect of the recent LHCb results for the 
decay B s — > \T on the Higgs sector, as well as the super-particle LHC searches and the dark 
matter constraints. Most importantly, we study other regimes than the decoupling regime of 
the pMSSM: the anti-decoupling regime for low Ma in which the roles of the h and H bosons 
are reversed, the intense coupling regime in which the three neutral particles h, A, H are rather 
close in mass, the intermediate regime at relatively low tan /3 in which the couplings of H, A 
to gauge+Higgs bosons are not too suppressed and the vanishing coupling regime in which 
the coupling of the h state to bottom quarks or gauge bosons is suppressed. Using the latest 
ATLAS, CMS and Tevatron data, we show that all these scenarios are now almost ruled out. 
Finally, we start studying the implications of the rates reported by the LHC experiments in the 
77 and ZZ final states used to obtain the Higgs boson signal and we comment on the bb final 
state to which the Tevatron is most sensitive. We perform a full scan of the pMSSM parameter 
space in order to delineate the regions which fit best the experimental data, including a possible 
enhancement of the h — > 77 rate. 

The paper is organised as follows. First, we briefly describe the pMSSM and its Higgs 
sector with its various regimes and summarise the Higgs decays and the production cross 
sections at the LHC. In section 3, we present the analysis of these different Higgs regimes 
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and the implications on the pMSSM parameters in the light of the LHC Higgs discovery and 
constraints. Section 4 has a short conclusion. 



2 The theoretical set-up 
2.1 The pMSSM Higgs sector 

In the MSSM the Higgs sector is extended to contain five Higgs particles Fj The lightest h 
boson has in general the properties of the Standard Model (SM) Higgs bosonjjand is expected 
to have a mass Mh ^ 115-135 GeV depending on the MSSM parameters, in particular, the 
ratio tan /3 of the vacuum expectation values of the two Higgs doublet fields that break the 
electro-weak symmetry in the MSSM. 

By virtue of supersymmetry, only two parameters are needed to describe the Higgs sector 
at tree-level. These can be conveniently chosen to be the pseudoscalar boson mass Ma and the 
ratio of vacuum expectation values of the two Higgs fields that break the symmetry, tan j3 = 
Vz/vi. However, accounting for the radiative corrections to the Higgs sector, known to play 



an extremely important role 1 10 , all soft SUSY-breaking parameters which are of (9(100) in 
addition to those of the SM, become relevant. This makes any phenomenological analysis in 
the most general MSSM a very complicated task. A phenomenologically more viable MSSM 
framework, the pMSSM, is defined by adopting the following assumptions: i) all soft SUSY- 
breaking parameters are real and there is no new source of CP-violation; ii) the matrices for 
the sfermion masses and for the trilinear couplings are all diagonal, implying no flavor change 
at tree-level; and ii) the soft SUSY-breaking masses and trilinear couplings of the first and 
second sfermion generations are the same at the electro- weak symmetry breaking scale. Making 
these three assumptions will lead to only 22 input parameters in the pMSSM: 

- tan/3: the ratio of the vevs is expected to lie in the range 1 < tan/3 < m t jm\;, 

- Ma- the pseudoscalar Higgs mass that ranges from M z to the SUSY-breaking scale; 

- fi: the Higgs-higgsino (supersymmetric) mass parameter (with both signs); 

- Mi, M2, M3: the bino, wino and gluino mass parameters; 

- mq,mu R ,rn^ R ,mi,me R \ the first/second generation sfermion mass parameters; 

- A u ,Ad,A e : the first/second generation trilinear couplings; 

- mq,mi R ,mi ,mi,m? R : the third generation sfermion mass parameters; 

- A t ,A b ,A T : the third generation trilinear couplings. 

Such a model has more predictability and it offers an adequate framework for phenomeno- 
logical studies. In general, only a small subset of the parameters appears when looking at a 
given sector of the pMSSM, such as the Higgs sector in this case. Some of these parameters will 
enter the radiative corrections to the Higgs boson masses and couplings. At the one-loop level, 
the h boson mass receives corrections that grow as the fourth power of the top quark mass 
m t (we use the running MS mass to re-sum some higher order corrections) and logarithmically 



1 For a review of the MSSM Higgs secto r [9] . For reviews on the radiative corrections in the MSSM Higgs 



sector and a complete set of references, see 10 
2 For a review of the SM Higgs boson, see 11 
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with the SUSY-breaking scale or common squark mass Ms', the trilinear coupling in the stop 



sector A t plays also an important role. The leading part of these corrections reads 12 



3 m\ 



2tt 2 v 2 sin 2 j3 



1 M l 
log^| + 
mf 



2 M 2 



XL 
6M| 



(1) 



We have defined the SUSY-breaking scale M s to be the geometric average of the two stop 
masses (that we take < 3 TeV not to introduce excessive fine-tuning) 

Ms = V™h^h ( 2 ) 
and introduced the mixing parameter X t in the stop sector (that we assume < 3Mg), 

X t = A t -fi,cotp. (3) 

The radiative corrections have a much larger impact and maximise the h boson mass in the 
so-called "maximal mixing" scenario, where the trilinear stop coupling in the DR scheme is 



maximal mixing scenario : 
In turn, the radiative corrections are much smaller for small values of X t , i.e. in the 



no mixing scenario : X t = 0. 



(4) 



(5) 



An intermediate scenario is when X+ is of the same order as Me which is sometimes called the 



typical mixing scenario : X t = Ms 



(6) 



These mixing scenarios have been very often used as benchmarks for the analysis of MSSM 
Higgs phenomenology 13 . The maximal mixing scenario has been particularly privileged since 
it gives a reasonable estimate of the upper bound on the h boson mass, M™ ax . We will discuss 
these scenarios but, compared to the work of Ref. fl~3], we choose here to vary the scale Ms- 



Together with the requirements on X t in eqs. (4-6), we adopt the following values for the 
parameters entering the pMSSM Higgs sector, 

1 



A t = A b , M 2 ~2M l 



-M s , M 3 = 0.8 M s . 



(7) 



and vary the basic inputs tan/3 and Ma- For the values tan/3 = 60 and Ma = Ms = 3 TeV 
and a top quark pole of mass of m t = 173 GeV, we would obtain a maximal Higgs mass value 
M™ ax Rj 135 GeV for maximal mixing once the full set of known radiative corrections up to 

In the no-mixing and typical mixing scenarios, one obtains 
D GeV and M™ ax « 125 GeV, respectively. Scanning over the 
soft SUSY-breaking parameters, one may increase these M™ ax values by up to a few GeV. 

It is important to note that the dominant two-loop corrections have been calculated in 
the DR scheme 



two loops is implemented 14 
much smaller values, M™ ax rj 



15 



16 



and SOFTSUSY 



19 



17 



that we 



and implemented in the codes Suspect 
will use here for the MSSM spectrum, but also in the on-shell scheme 18 as implemented in 
FeynHiggs 



In general, the results for Mh in the two scheme differ by at most 2 GeV, 
which we take as a measure of the missing higher order effects. Quite recently, the dominant 
three-loop contribution to Mh has been calculated and found to be below 1 GeV [20]. Thus, 
the mass of the lightest h boson can be predicted with an accuracy of AMh ~ 3 GeV and this 
is the theoretical uncertainty on Mh that we assume. 
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2.2 The various regimes of the pMSSM 

The spectrum in the various regimes of the pMSSM Higgs sector M, depends on the values of 
Ma and also on tan/3, and that we will confront to the latest LHC and Tevatron data in this 
paper. 



We start from the decoupling regime 21 that has been already mentioned and which in 
principle occurs for large values of M4 but is reached in practice at Ma ?t 300 GeV for low 
tan/3 values and already at Ma ^ M™ ax for tan/3 > 10. In this case, the CP-even h boson 
reaches its maximal mass value M™ ax and its couplings to fermions and gauge bosons (as well 
as its self-coupling) become SM-like. The heavier H boson has approximately the same mass 
as the A boson and its interactions are similar, i.e. its couplings to gauge bosons almost vanish 
and the couplings to bottom (top) quarks and r leptons fermions are (inversely) proportional 
to tan/3. Hence, one will have a SM-like Higgs boson h = Hsm and two pseudo-scalar (like) 
Higgs particles, $ = H, A. The H boson is also degenerate in mass with the A boson and 
the intensity of its couplings to fermions is similar. Hence, in the decoupling limit, the heavier 
H/A/H ± bosons almost decouple and the MSSM Higgs sector reduces effectively to the SM 
Higgs sector, but with a light h boson. 



The anti-decoupling regime 22 occurs for a light pseudo-scalar Higgs boson, Ma ^ M™ a> 



L h i 

and is exactly opposite to the decoupling regime. The roles of the h and H bosons are reversed 

and at large tan (3 values, the h boson is degenerate in mass with the pseudo-scalar A, Mh — 

Ma, while the H boson has a mass close to its minimum which is in fact M™ ax . Because of the 

upper bound on M h , all Higgs particles are light. Here, it is the h boson which has couplings 

close to those of A, $ = h, A, while the H boson couplings are SM-like, H = H SM . 



The intense-coupling regime 23 occurs when the mass of the pseudo-scalar A boson is close 
to M™ ax . In this case, the three neutral Higgs bosons h, H and A (as well as the charged Higgs 
particles) have comparable masses, Mh ~ Mh ~ Ma ~ M™ ax . The mass degeneracy is more 
effective when tan /3 is large. Here, both the h and H bosons have still enhanced couplings 
to 6-quarks and r leptons and suppressed couplings to gauge bosons and top quarks, as is 
the pseudo-scalar A. Hence, one approximately has three pseudo-scalar like Higgs particles, 
$ = h, H, A with mass differences of the order of 10-20 GeV. 

The intermediate-coupling regime occurs for low values of tan /3, tan f3 < 5-10, and a not 



too heavy pseudo-scalar Higgs boson, Ma ^ 300-500 GeV [9|. Hence, we are not yet in the 
decoupling regime and both CP-even Higgs bosons have non-zero couplings to gauge bosons 
and their couplings to down-type (up-type) fermions (as is the case for the pseudoscalar A 
boson) are not strongly enhanced (suppressed) since tan /3 is not too large. This scenario is 
already challenged by LEP2 data which call for moderately large values of tan /3. 

The vanishing-coupling regime occurs for relatively large values of tan /3 and intermediate 
to large Ma values, as well as for specific values of the other MSSM parameters. The latter 
parameters, when entering the radiative corrections, could lead to a strong suppression of the 
couplings of one of the CP-even Higgs bosons to fermions or gauge bosons, as a result of the 



cancellation between tree-level terms and radiative corrections 24 . An example of such a 
situation is the small a e fr scenario which has been used as a benchmark [13] and in which the 
Higgs to bb coupling is strongly suppressed. 

Within the plane [Ma, tan /3] , the parameter space in which the above regimes of the pMSSM 
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Figure 1: The parameter space for the various regimes of the MSSM Higgs sector as defined in the 
text and in eq. (8) in the [Ma, tan/3] plane; the maximal mixing scenario with Ms = 2 TeV is adopted. 



Higgs sector occur are displayed in Figure [TJ We have chosen the usual maximal mixing scenario 
with Ms = 2 TeV and the other SUSY parameters as in eq. ([7]), except for the vanishing 
coupling scenario, where we have scanned over the SUSY parameters, and only w 5 x 10~ 4 of 
the scanned points fulfil its requirements. The following conditions have been imposed: 

cos 2 (/3-a) < 0.05 
cos 2 (/3-a) > 0.95 

0.05 < cos 2 (/3 - a) < 0.7, tan/3 < 10 
M A < 140 GeV, g 2 hbb and g 2 Hbb > 50 
M A > 200 GeV, g 2 hbb or g 2 hVV < 0.05. (8) 
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In addition, we have to consider the SUSY regime, in which some SUSY particles such as 
the charginos, neutralinos as well as the third generation sleptons and squarks, could be light 
enough to significantly affect the phenomenology of the MSSM Higgs bosons. For instance, 
light sparticles could substantially contribute to the loop induced production and decays modes 



of the lighter h boson 25 , 26 and could even appear (in the case of the lightest neutralino) in 



its decay product as will be discussed below. 



2.3 Higgs decays and production in the pMSSM 

For the relatively large values of tan/3 presently probed at the LHC, tan/3 > 7 as discussed 
below, the couplings of the non-SM like Higgs bosons to b quarks and r leptons are so strongly 
enhanced and those to top quarks and gauge bosons suppressed, that the pattern becomes as 
simple as the following (more details can be found in Ref. [9]): 

- The $ = A or H/h bosons in the decoupling/anti-decoupling limit decay almost exclu- 
sively into bb and t + t~ pairs, with branching ratios of, respectively, ~ 90% and ~ 10%, and 
all other channels are suppressed to a level where their branching ratios are negligible. 
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- The H ± particles decay into fermion pairs: mainly H + — > tb and H + — > ru T final states 
for H ± masses, respectively, above and below the tb threshold. 

- The CP-even h or H boson, depending on whether we are in the decoupling or anti- 
decoupling regime, will have the same decays as the SM Higgs boson. For M h / H « 126 GeV, 
the main decay mode will be the bb channel with a ~ 60% probability, followed by the decays 
into cc, t + t~ and the loop induced decay into gluons with ~ 5% branching ratios. The WW* 
decay reaches the level of 20%, while the rate for ZZ* is a few times 10~ 2 . The important loop 
induced 77 decay mode which leads to clear signals at the LHC have rates of 0(1O~ 3 ). 

In the intense-coupling regime, the couplings of both h and H to gauge bosons and up-type 
fermions are suppressed and those to down-type fermions are enhanced. The branching ratios 
of the h and H bosons to bb and t + t~ final states are thus the dominant ones, with values 
as in the case of the pseudoscalar A boson. In the intermediate-coupling regime, interesting 
decays of H, A and H into gauge and/or Higgs bosons occur, as well as A/H — )■ ti decays, but 
they are suppressed in general. Finally, for the rare vanishing-coupling regime when the Higgs 
couplings to 6-quarks and eventually r-leptons accidentally vanish, the outcome is spectacular 
for the h boson: the WW* mode becomes dominant and followed by h — > gg, while the 
interesting h — > 77 and h — > ZZ* decay modes are enhanced. 

In the case of the SM-like Higgs particle (that we assume now to be the h boson), there 
are two interesting scenarios which might make its decays rather different. First we have the 
scenario with the Higgs bosons decaying into supersymmetric particles. Because most sparticles 
must be heavier than about 100 GeV, there is no SUSY decays of the h boson except for the 
invisible channel into a pair of the lightest neutralinos, h — > XxXi- This is particularly true 
when the gaugino mass universality relation M 2 ~ 2Mi is relaxed, leading to light x± states 
while the LEP2 bound, m v ± > 100 GeV, still holds. In the decoupling limit, the branching 
ratio of the invisible decay can reach the level of a few 10%. Decays of the heavier Aj H/ H 
bosons, in particular into charginos, neutralinos, sleptons and top squarks, are in turn possible. 
However, for tan/3 > 10, they are strongly suppressed. 

The second scenario of interest occurs when SUSY particles contribute to loop-induced 
Higgs decays. If scalar quarks are relatively light, they can lead to sizable contributions to the 
decays h — >• gg and h — > 77. Since scalar quarks have Higgs couplings that are not proportional 
to their masses, their contributions are damped by loop factors l/ m /> and decouple from the 
vertices contrary to SM quarks. Only when nig is not too large compared to that the 



contributions are significant 25 . This is particularly true for the ti contributions to h — > gg, 
the reasons being that large X t mixing leads to a t\ that is much lighter than all other squarks 
and that the h coupling to stops involves a component which is proportional to m t X t and, for 
large X t , it can be strongly enhanced. Sbottom mixing, oc m^X^, can also be sizable for large 
tan/3 and /i values and can lead to light 61 states with strong couplings to the h boson. In 
h — > 77 decay, there are in addition slepton loops, in particular f states which behave like 
scalar bottom quarks and have a strong mixing at high \i tan (3, can make a large impact on the 
decay rate. Besides, chargino loops also enter the h — > 77 decay mode but their contribution 
is in general smaller since the Higgs-xx couplings cannot be strongly enhanced. 

For the evaluation of the decay branching ratios of the MSSM Higgs bosons, we use the 



program HDECAY 29 , which incorporates all decay channels including those involving super- 



particles and the most important sets of higher order corrections and effects. 
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Coming to Higgs boson production at the LHC, for a SM-like particle H^m there are 



essentially four mechanisms for single production 11 . These are gg fusion, gg — > H^m, vector 



boson fusion, qq — > HsmQQ, Higgs-strahlung, qq — > HsmV and tt associated Higgs production, 
pp — > tiHsM- The gg — > H$m process proceeds mainly through a heavy top quark loop and 
is by far the dominant production mechanism at the LHC. For a Higgs boson with a mass of 
~ 126 GeV, the cross section is more than one order of magnitude larger than in the other 
processes. Again for M# SM ~ 126 GeV, the most efficient detection channels are the clean but 
rare H — >■ 77 final states, the modes H — > ZZ* — > 4£ ± , H — > WW^ — > ttvv with t = e,/i 
and, to a lesser extent, also Hsm — > t + t~ . At the LHC and, most importantly, at the Tevatron 
one is also sensitive to qq — > H$m + W/Z — > bb + W/Z with W — )■ tv and Z — > ££, vv. 

For the MSSM Higgs bosons, the above situation holds for the h(H) state in the (anti-) 
decoupling regime. Since AVV couplings are absent, the A boson cannot be produced in 
Higgs-strahlung and vector boson fusion and the rate for pp — > tiA is strongly suppressed. 
For the $ = A and h(H) states, when we are in the (anti-) decoupling limit, the b quark will 
play an important role for large tan/3 values as the $66 couplings are enhanced. One then 
has to take into account the 6-loop contribution in the gg — >• $ processes which becomes the 
dominant component in the MSSM and consider associated Higgs production with bb final 
states, pp — > 66 + $ which become the dominant channel in the MSSM. The latter process is in 
fact equivalent to 66 — > $ where the 6-quarks are taken from the proton in a five active flavor 
scheme. As the $ bosons decay mainly into 66 and t + t~ pairs, with the former being swamped 
by the QCD background, the most efficient detection channel would be pp — > $ — > t + t~ . This 
process receives contributions from both the gg — >• $ and 66 — > $ channels. 

These processes also dominate the h/H/A production in the intense coupling regime. In 
fact, in the three regimes above, when all processes leading to r + r~ final states are added up, 
the rate is 2 x a(gg + 66 — >■ A) x BR(A — > r + r _ ). In the intermediate coupling regime, these 
process have very low cross sections as for 3-5 < tan/3 < 7-10, the $66 couplings are not 
enough enhanced and the $££ ones that control the gg fusion rate are still suppressed. 

Finally, for the charged Higgs boson, the dominant channel is the production from top 
quark decays, t — > H + b, for masses not too close to M H ± = m t —mb. This is true in particular 
at low or large tan /3 values when the t — > H + b branching ratio is significant. 

The previous discussion on MSSM Higgs production and detection at the LHC might be 
significantly altered if scalar quarks, in particular t and 6, are light enough. Indeed, the Hgg 
and hgg vertices in the MSSM are mediated not only by the t/b loops but also by loops 
involving their partners similarly to the Higgs photonic decays. The gg — )■ h cross section in 
the decoupling regime can be significantly altered by light stops and a strong mixing X t which 
enhances the htiti coupling. The cross section times branching ratio u{gg — > h) x BR(/i — > 77) 
for the lighter h boson at the LHC could be thus different from the SM case, even in the 



decoupling limit in which the h boson is supposed to be SM-like 25 



Finally, we should note that in the scenario in which the Higgs bosons, and in particular 
the lightest one h, decay into invisible lightest neutralinos, h — > X1X11 the observation of the 
final state will be challenging but possible at the LHC with a higher energy and more statistics. 
This scenario has recently been discussed in detail in Refs. [5 30 . 



S 



3 Analysis and results 



3.1 pMSSM scans and software tools 

The analysis is based on scans of the multi-parameter MSSM phase space. The input values 
of the electro-weak parameters, i.e. the top quark pole mass, the MS bottom quark mass, the 
electro-weak gauge boson masses, electromagnetic and strong coupling constants defined at 
the scale Mz, are given below with their la allowed ranges |32|, 

m t = (173 ± 1) GeV, m b (m b ) = (4.19^) GeV, 
M z = (91.19 ± 0.002) GeV, M w = (80.42 ± 0.003) GeV, 
a(M|) = 1/127.916 ± 0.015, ot s (M 2 z ) = 0.1184 ± 0.0014 . (9) 

The pMSSM parameters are varied in an uncorrelated way in flat scans, within the following 
ranges: 

1 < tan/3 < 60, 
50 GeV < M A < 3 TeV , 
-10 TeV < A f < 10 TeV, 
50 GeV < m f - L ,m h ,M 3 < 3.5 TeV, 

50 GeV < Mi, M 2 , < 2.5 TeV (10) 

to generate a total of 6 x 10 7 pMSSM points. The scan range is explicitly chosen to include 
the various mixing scenarios in the Higgs section discussed in section 2.1: the maximal mixing, 
no-mixing and typical mixing scenarios. Additional 10 7 points are generated in specialised 
scans used for the studies discussed later in section 3.5. We select the set of points fulfilling 
constraints from flavour physics and lower energy searches at LEP2 and the Tevatron, as 
discussed in Ref. [33], to which we refer also for details on the scans. We highlight here the 



tools most relevant to this study. The SUSY mass spectra are generated with SuSpect 16 



and SOFTSUSY 3.2.3 [17] . The superparticle partial decay widths and branching fractions are 



computed using the program SDECAY 1.3 34 . The flavour observables and dark matter relic 



density are calculated with Super I so Relic v3.2 35. 

The Higgs production cross sections at the LHC are computed using HIGLU 1.2 |36| for 
the gg — >• h/H/A process, including the exact contributions of the top and bottom quark loops 
at NLO-QCD and the squark loops, and the program bbOnnlo for bb — > h/H/A at NNLO- 
QCD. They are interfaced with Suspect for the MSSM spectrum and HDECAY for the Higgs 
decay branching ratios. The Higgs production cross sections and the branching fractions for 
decays into bb, 77, WW and ZZ from HIGLU and HDECAY are compared to those predicted by 
FeynHiggs. In the SM both the gg — > H$m cross section and the branching fractions agree 
within ~ 3%. Significant differences are observed in the SUSY case, with HDECAY giving values 
of the branching fractions to 77 and WW, ZZ which are on average 9% lower and 19% larger 
than those of FeynHiggs and have an r.m.s. spread of the distribution of the relative difference 
between the two programs of 18% and 24%, respectively (H]. 
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3.2 Constraints 



We apply constraints from flavour physics, anomalous muon magnetic moment, dark matter 
constraints and SUSY searches at LEP and the Tevatron. These have been discussed in details 
in Ref. [33]. In particular, we consider the decay B s — > which can receive extremely large 

SUSY contributions at large tan/3. An excess of events in this channels has been reported by 
the CDF-II collaboration at the Tevatron [37] and upper limits by the LHCb [38] and CMS [39] 
collaborations at LHC Recently the LHCb collaboration has presented their latest result for 
the search of this decay based on 1 fb -1 of data. A 95% C.L. upper limit on its branching 



fraction is set at 4.5 x 10 9 38 . After accounting for theoretical uncertainties, estimated at 
the 11% level |40] the constraint 



BR(B S n+fT) < 5 x 10~ 9 (11) 
is used in this analysis. For large values of tan/3, this decay can be enhanced by several orders 



of magnitude so that strong constraints on the scalar contributions can be derived 41 , and the 



small Ma and large tan (3 region can be severely constrained. As already remarked in Ref. 33 



the constraints obtained are similar and complementary to those from the dark matter direct 



detection limits of XENON-100 42 and searches for the A — > t + t decay. 



Concerning the relic density constraint, we impose the upper limit derived from the WMAP- 



7 result 43 



10" 4 < tt x h 2 < 0.155 , (12) 



accounting for theoretical and cosmological uncertainties (44] . 

The searches conducted by the ATLAS and CMS collaborations on the y/s = 7 TeV data 



for channels with missing Et 4546 have already provided a number of constraints relevant 
to this study. These have excluded a fraction of the pMSSM phase space corresponding to 
gluinos below ~ 600 GeV and scalar quarks of the first two generations below ~ 400 GeV. 



These constraints are included using the same analysis discussed in Ref. 33 , extended to an 
integrated luminosity of 4.6 fb _1 . 



Then, searches for the MSSM Higgs bosons in the channels h/H/A — > t + t~ 47 , 48 have 
already excluded a significant fraction of the [Ma, tan /3] plane at low Ma values, Ma ^ 200 
GeV and tan /3 < 10, and larger values of tan f3 for Ma ^ 200 GeV. These constraint on the 
pMSSM parameter space are already important. It is supplemented by the search of light 
charged Higgs bosons in top decays, t — > bH + — > bru, performed by the ATLAS collaboration 



49 which is effective at low Ma values, Ma ^ 140 GeV, corresponding to M H ± < 160 GeV. 

Following the Higgs discovery at the LHC, the lightest Higgs boson in our analysis is 
restricted to have a mass in the range allowed by the results reported by ATLAS and CMS: 

123 GeV < M h < 129 GeV (13) 

where the range is centred around the value corresponding to the average of the Higgs mass 
values reported by ATLAS and CMS, Mh — 126 GeV, with the lower and upper limits ac- 
counting for the parametric uncertainties from the SM inputs given in eq. ([9]), in particular the 
top quark mass, and the theoretical uncertainties in the determination of the h boson mass. 
It is also consistent with the experimental exclusion bounds. 
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The impact of the Higgs mass value and its decay rates on the parameters of the pM SSM can 
be estimated by studying the compatibility of the pMSSM points with the first results reported 
by ATLAS [I] and CMS [2] at the LHC and also by the Tevatron experiments [3] . Starting from 
our set of 6 x 10 7 pMSSM points which are pre-selected for compatibility with the constraints 
discussed above, we consider the two decay channels giving the Higgs boson evidence at the 
LHC, 77 and ZZ and include also the bb and rr channels. In the following, we use the notation 
Rxx to indicate the Higgs decay branching fraction to the final state XX, BR(/i — > XX), 
normalised to its SM value. We also compute the ratios of the product of production cross 
sections times branching ratios for the pMSSM points to the SM values, denoted by \ixx 
for a given h — > XX final state, fixx = ff (^jxBR(%H^) ' These are com P are d to the 
experimental values. For the 77, and ZZ channels we take a weighted average of the results 
just reported by the experiments, as given in Table [T] with their estimated uncertainties. 



Parameter 


Value 


Experiment 


M H 


125.9±2.1 GeV 


ATLAS 


I 


+ CMS 


IL 




1.71±0.33 


ATLAS 


50 


+ CMS 


13 


Hzz 


0.95±0.40 


ATLAS 


52 


+ CMS 






<1.64 (95% C.L.) 


CMS m 






fJ>rr 


<1.06 (95% C.L.) 


CMS |55 





Table 1: Input parameters used for the pMSSM study. 



While the results are compatible with the SM expectations within the present accuracy, 
they highlight a possible enhancement in the observed rates for the 77 channel, where ATLAS 
and CMS obtain /i 77 = 1.9±0.5 and 1.56±0.43, respectively. In the following, we do not take 
into account the theoretical uncertainties in the production cross section, which are estimated 



significant for the main production channel, gg — >• h 56 57 



3.3 The decoupling regime 

Figure [2] presents the parameter space [Ma, tan /5] in our benchmark scenario with M5 = 2 TeV 
in the maximal mixing scenario. The regions excluded by the various constraints that we have 
imposed are indicated. The green area corresponds to the non-observation of Higgs bosons at 
LEP2 which excludes tan/3 < 3 at moderate to large Ma values, Ma ^ 150 GeV, but up to 
tan /5 rj 5-10 at low Ma values. The blue area is the one ruled out by the latest published 
results of the CMS collaboration on the search of resonances decaying into r + r~ final states; 
it touches the LEP2 band at small Ma, but reduces in size when Ma is increased. The small 
visible area in red is the one excluded by the B s — > constraint but, in fact, part of the 

excluded region is hidden by the CMS blue area. 

To that, we superimpose the area in which we make the requirement 123 < Mh < 129 GeV, 
that is indicated in dark blue. This band covers the entire range of Ma values and leaves only 
the tan/5 values that are comprised between tan/5 ~ 3-5 and tan/5 ~ 10. Between the 
LEP2 and the "M^" blue band, one has Mh < 123 GeV, while above the Mh band, one has 
Mh > 129 GeV and both areas are excluded. The requirement that the h boson mass should 
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have the value measured at the LHC, even with the large uncertainty that we assume, provides 
thus a strong constraint on the [Ma, tan/3] parameter space in the pMSSM. 



Maximal mixing • M g = 2 TeV 




100 150 200 250 300 350 400 450 500 

M A (GeV) 

Figure 2: The parameter space [Ma, tan f3] for Ms = 2 TeV in the maximal mixing scenario with the 
individual constraints from LEP2 (green), CMS t + t~ searches (light blue) and flavor physics (red) 
displayed. The area in which 123 < Mh < 129 GeV is also shown (dark blue). 



In Figure [3j we show the same [Ma, tan /3] plane but for different SUSY-breaking scales, 
Ms = 1,2 and 3 TeV and for the zero, typical and maximal mixing scenarios defined in eqs. (4- 
6). As can be seen, the situation changes dramatically depending on the chosen scenario. Still, 
in the maximal mixing scenario with Ms = 3 TeV the size of the Mh band is reduced from 
above, as in this case, already values tan /3 > 5 leads to a too heavy h boson, Mh > 129 GeV. 
In turn, for Ms = 1 TeV, the entire space left by the LEP2 and CMS Higgs constraints is 
covered with many points at tan/3 > 20 excluded by the flavor constraint. Nevertheless, the 
possibility with Ms ~ 1 TeV will start to be challenged by the search for squarks at the 
LHC when 30 fb _1 of data will be collected by the experiments. In the no-mixing scenario, 
it is extremely hard to obtain a Higgs mass of Mh > 123 GeV and all parameters need to be 
maximised: Ms = 3 TeV and tan /3 > 20; a small triangle is thus left over, the top of which is 
challenged by the flavor constraints. The typical mixing scenario resembles to the no-mixing 
scenario, with the notable difference that for Ms = 3 TeV, the entire space not excluded by 
the LEP2 and CMS constraints allow for an acceptable value of Mh- 

In the discussion so far, we have adopted the value m t = (173±1) GeV for the top quark 
mass as measured by the CDF and DO experiments at the Tevatron 58 . This implicitly 



assumes that this mass corresponds to the top quark pole mass, i.e. the mass in the on-shell 
scheme, which serves as input in the calculation of the radiative corrections in the pMSSM 
Higgs sector and, in particular, to the mass Mh- However, the mass measured at the Tevatron 
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Maximal mixing - M = 1 TeV 




Maximal mixing - M = 2 TeV 



Maximal mixing - M = 3 TeV 
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Typical mixing - M g = 1 TeV 
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Figure 3: The [M^,tan/3] plane for M$ = 1, 2 and 3 TeV and for zero, typical and maximal mixing. 
The colour coding for the different regions is the same as in Fig. [2] 



is not theoretically well denned and it is not proved that it corresponds indeed to the pole mass 
as discussed in 59 . For an unambiguous and well-defined determination of the top quark mass, 



it is appears to be safer to use the value obtained from the determination of the top quark pair 
production cross section measured at the Tevatron, by comparing the measured value with 
the theoretical prediction at higher orders. This determination has been recently performed 
yielding the value of (173. 3 ±2. 8) GeV 59 for m^ ole . The central value is very close to that 
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measured from the event kinematics but its uncertainty is larger as a result of the experimental 
and theoretical uncertainties that affect the measurement. 

It is interesting to assess the impact of a broader mass range for the top quark. We return 
to our benchmark scenario with M s = 2 TeV and maximal stop mixing and draw the "M^" 
bands using the top quark mass values of 170 GeV and 176 GeV corresponding to the wider 
uncertainty interval quoted above. The result is shown in Figure |4j A 1 GeV change in m t input 
leads to a ~1 GeV change in the corresponding Mh value. The smaller value of m t would open 
up more parameter space as the region in which Mh > 129 GeV will be significantly reduced. 
In turn, for m t = 176 GeV, the corresponding h boson mass increases and the dark-blue 
area quite significantly shrinks, as a result. It must be noted that for m t = 170 GeV, the 
no-mixing scenario would be totally excluded for Ms < 3 TeV, while in the typical mixing 
scenario only a small area at high tan /3 will remain viable. For m t = 176 GeV significant 
[Ma, tan fl] regions that was excluded when taking the ± 1 GeV uncertainty for top mass 
value becomes allowed. The impact of the value of mt is thus extremely significant. This is 




- -I h I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i H I- 1 

100 150 200 250 300 350 400 450 500 100 150 200 250 300 350 400 450 500 100 150 200 250 300 350 400 450 500 

M A (GeV) M A (GeV) M A (GeV) 



Figure 4: The parameter space [Ma, tan 0\ for Ms = 2 TeV, maximal mixing and three values of the 
top quark mass mt = 170 GeV (left), 173 GeV (centre) and 176 GeV (right). 

even more true in constrained scenarios, where the top mass also enters in the evaluation of 
the soft SUSY-breaking parameters and the minimisation of the scalar potential. To visualise 
the impact of m t , we have repeated the study presented in Ref. j2], presenting the maximal 
Mh value reached when scanning over all the parameters of the minimal SUGRA, AMSB and 
GMSB models. Figure [5] shows the result with the M™ ax value as a function of Ms taking 
m t =173±3 GeV. While for m t = 173 GeV, there is no region of the parameter space of the 
mAMSB and mGMSB models which satisfies M h > 123 GeV, for M s < 3 TeV assumed in g], 
and the models are disfavoured, using m t =176 GeV, the regions of these mAMSB and mGMSB 
models beyond Ms = 2 TeV become again viable. This will be also the case of some of the 
variants and even more constrained mSUGRA scenarios. Further, even for m(=173 GeV, if we 
move the Ms upper limit from the 3 TeV boundary adopted in Ref. R] to Ms = 5 TeV, these 
models have region of their parameters compatible with the LHC Higgs mass. 

Finally, we comment on the impact of increasing the Mh allowed range from 123 GeV 
< Mh < 127 GeV as was done in Ref. [4] relying on the 2011 LHC data, to the one adopted here, 
123 GeV < Mh < 129 GeV, in the various constrained models discussed in that reference (and 
to which we refer for the definition of the models and for the ranges of input parameters that 
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Figure 5: Maximal Higgs mass in the constrained MSSM scenarios mSUGRA, mAMSB and mGMSB, 
an a function of the scale Mg when the top quark mass is varied in the range nit = 170-176 GeV. 



have been adopted). The outcome is shown in Fig. |6j where the maximal h mass value obtained 
by scanning the basic input parameters of the model over the appropriate ranges. In the left- 
hand side, M™ ax is displayed as a function of tan and in the right-hand side as a function 
of Mg. As the lower bound M™ ax > 123 GeV is the same as in our previous analysis, the 
mASMB, mGMSB and some variants of the mSUGRA model such as the constrained NMSSM 
(cNMSSM), the no-scale model and the very constrained MSSM (VCMSSM) scenarios are still 
disfavoured. However, for mSUGRA and the non-universal Higgs mass model (NUHM), all 
values of tan/3 > 3 and 1 TeV < M s < 3 TeV lead to an appropriate value of M h when 
including the uncertainty band. 



•135 
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Figure 6: The maximal h mass value M™ ax as functions of tan /3 (left) and Mg (right) in the mASMB, 
mGMSB as well as in mSUGRA and some of its variants. The basic parameters of the models are 
varied within the ranges given in Ref. (4j; the top quark mass is fixed to mt = 173 GeV. 
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3.4 The other (non— SUSY) regimes 



The other regimes of the pMSSM Higgs sector, apart from the decoupling and the SUSY 
regimes, occur for low to intermediate values of the pseudoscalar Higgs mass, Ma ^ 200 GeV, 
and relatively large tan values. These are the anti-decoupling, the intense, the intermediate 
and the vanishing coupling regimes. These are constrained by the results of the LEP2 and 
LHC searches. The LEP2 results for Ma ^ 200 GeV and not too large M s values, lead to 
tan (3 > 3, 8 and 10 for, respectively the maximal, typical and no-mixing scenarios; see Fig. |3j 
The negative search for Higgs particles in r-lepton final states, pp — > $ — > t + t~, by the 
ATLAS and CMS collaborations places further constraints. While in the decoupling regime, 
the relevant Higgs states would be $ = A + H, these are $ = A + h and $ = A + H + h 
in the anti-decoupling and intense coupling regimes, respectively. As already mentioned, one 
would have in the three regimes the same signal cross section times branching ratios u{pp — >• 
$ — > r + r~) f=s a(bb + gg — >■ A) x BR(A — > t + t~) almost independently of the mixing scenario 
and the other pMSSM parameters [57] . The constraint from the CMS published results alone 
with the ~ 5 fb _1 of data collected in 2011 48 imposes tan/3 > 10, as shown in Fig. [7] where 
we zoom on the plane [Ma, tan /3] at low to intermediate Ma values, for the maximal mixing 
scenario and Ms = 2 TeV. 

This limit can be strengthened by the same t + t~ search performed by the ATLAS col- 
and also by the t — > H + b search in top decays [49 which is effective for 
130 GeV and which, as can be seen in Fig. 7| excludes large tan/3 



47 



H± 



+ Ml < 



laboration 

M A 

values for which BR(t — > bH + ) is significant. Put together, these constraints exclude entirely 
both the anti-decoupling and intense coupling regimes. Would remain then, the intermediate 
coupling regime with tan /3 5-8 when the LEP2 constraint is also imposed. Depending on 
the mixing scenario, most of it will be excluded by the Mh ~ 126 GeV constraint (see Fig. |3|). 



M s = 2 TeV - Maximal mixing 
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I I intermediate 
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ntense 
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Figure 7: Parameter space for the various regimes of the MSSM Higgs sector as defined in the text 
and in eq. (8) in the [M4,tan/3] plane, in the maximal mixing scenario with Ms = 2 TeV. The upper 
limit constraints from $ — > tt (continuous light blue line) and t — > H + b (dashed blue line) searches 
at the LHC are shown together with the LEP2 excluded region (green area). 

A very interesting possibility would be that the observed Higgs particle at the LHC is 
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actually the H state, while the lighter h boson has suppressed couplings to W/Z bosons and 
top quarks, allowing it to escape detection. In this case, the H couplings to bottom quarks 
should not be too enhanced, tan/3 < 8, not to be in conflict with the r + r~ and t — > bH + 
searches above. For the H boson to be SM-like, one should have Ma ~ Mh ~ 126 GeV 
and not too low tan/3 values, tan /3 > 7-10. One is then in the borderline between the anti- 
decoupling and the intermediate coupling regimes. We have searched for points in which indeed 
M H ~ 126 GeV with couplings to VV states, gnvv ~ 0.9, such that the H — > ZZ and H — > 77 
(which mainly occurs through a H^-boson loop) decays are not suppressed compared to the 
measured values by ATLAS and CMS given in Table HI In our scan, out of the 10 6 points, 
before imposing any LHC-Higgs constraint, only « 20 points fulfilled the above requirements. 
These points are then completely excluded once the flavor constraints, in particular those from 
the b — > 57 radiative decay, are imposed. Hence, the possibility that the observed Higgs particle 
at the LHC is not the lightest h particle appears highly unlikely according to the result of our 
scan of the parameter space. Combining the h/A — >■ t + t~ and the t — > bH + constraints and 
including the results on the new 8 TeV data should further constrain the parameter space and 
completely exclude this scenario. 

Finally, the vanishing coupling regime is strongly disfavoured by the LHC and Tevatron 
data that are summarised in Table [T} The observation of H — > ZZ final states by the ATLAS 
and CMS collaborations rules out the possibility of vanishing hVV couplings. The reported 
excess of events in the qq — > VH — > Vbb process by the CDF and DO collaborations seem also 
to rule out both the vanishing hbb and hVV coupling scenarios. However, there is still the 
possibility that these couplings are smaller than those predicted in the SM case, in particular 



because of the effects of SUSY particles at high tan /3 24 . We are then in the SUSY-regime 
to which we turn now. 



3.5 The SUSY regime 

In the SUSY regime, both the Higgs production cross section in gluon-gluon fusion and the 
Higgs decay rates can be affected by the contributions of SUSY particles. This makes a 
detailed study of the pMSSM parameter space in relation to the first results reported by the 
ATLAS and CMS collaborations especially interesting for its sensitivity to specific regions of 
the pMSSM parameter space. In particular, the branching fraction for the 77 decay of the h 
state is modified by Higgs mixing effects outside the decoupling regime as was discussed above, 
by a change of the hbb coupling due to SUSY loops [24] , by light superparticle contributions 
to the /177 vertex (6 25 26 and by invisible h decays into light neutralinos |28|. 



We study these effects on the points of our pMSSM scan imposing the LHC results as 
constraints. The numerical values adopted in the analysis are given in Table [TJ assuming in 
the following on that the observed particle is the h state. First, we briefly summarise the 
impact of the SUSY particles on the Higgs decay branching fractions, staring from invisible 
decays, and production cross sections. Then we discuss our finding on the impact of the LHC 
and Tevatron data on the pMSSM parameters. 
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3.5.1 Invisible Higgs decays 



Despite the fact that the discovered particle has a sufficient event rate in visible channels to 
achieve its observation, it is interesting to consider the regions of parameter space in which 
invisible Higgs decays occur. This scenario has recently been re-considered in (5 30 . Besides 
the value of M^, the invisible branching ratio BR(/i — > XiXi) is controlled by four parameters: 
the gaugino masses Mi and M 2 , the higgsino parameter \i and tan/3. They enter the 4x4 
matrix Z which diagonalises the neutralino mass matrix. They also enter the Higgs coupling 
to neutralinos which, in the case of the LSP, is 



9hx lX iO ^ ( z i2 - tan%Zn) (sin (3Z 13 + cos (3Z 



u 



(14) 



if we assume the decoupling limit not to enhance the h — )■ 66 channel which would significantly 
reduce the invisible decay. In this coupling, Z n , Z 12 are the gaugino components and Z ri , Z u 
the higgsino components. Thus, the coupling vanishes if the LSP is a pure gaugino, ^> Mi 
leading to m x o « Mi, or a pure higgsino, Mi S> |//| with m x o ~ 

For the invisible decay to occur, a light LSP, m x o < \Mh is required. Since in the pMSSM, 
the gaugino mass universality M 2 ~ 2Mi is relaxed, one can thus have a light neutralino 
without being in conflict with data. The constraint from the Z invisible decay width measured 
at LEP restricts the parameter space to points where the Xi is bino-like, if its mass is below 
45 GeV, and thus to relatively large values of the higgsino mass parameter Since a large 
decay width into XiXi corresponds to small values of this remove a large part of the 
parameter space where the invisible Higgs decay width is sizable. Still, we observe invisible 
decays for 45 GeV< M x o < M h o/2 and < 150, corresponding to a combination of parameters 
where the x? is a mixed higgsino-gaugino state |5|. These pMSSM points are shown in the 
[Mi, /i] plane in the left panel of Fig. M 
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Figure 8: Left: Points in the [Mi, fj] parameter space where the invisible branching fraction BR(/i — > 
XiXi) ^ 0.15 from a pMSSM scan where we impose the LEP constraint on the Z invisible width and 
neutralino relic density Q c hih 2 . Right: r2 x /i 2 as a function of m^o with all the selected pMSSM points 
in black and those giving a BR(/i — > XiXi) ^ 0-15 in grey. 



If the LSP at such a low mass were to be the dark-matter particle, with the relic density 
given in eq. (12), it should have an efficient annihilation rate into SM particles. The only 
possible way for that to occur would be XiXi annihilation through the s-channel light h 
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pokj^ 60 which implies that m x o < to still have a non-zero invisible branching ratio, as 



shown in the right panel of Fig. 8j where the pMSSM points satisfying BR(/i — » > 5% 

are shown in the plane [m x oAog 10 (Qh 2 )}. However, because the partial decay width T(h — >■ 
XiXi) is suppressed by a factor /3 3 near the Mh ~ 2m x o threshold, with the velocity /3 = 
(1 — 4m 2 /M|) 1 / 2 , the invisible branching fraction is rather small if the WMAP dark matter 
constraint is to hold. MSSM light neutralinos compatible with claims of direct detection dark 
matter signals are also consistent with collider bounds |61|. 



3.5.2 Sparticle effects on the hbb coupling 

SUSY particles will contribute to the hbb coupling as there are additional one-loop vertex cor- 



rections that modify the tree-level Lagrangian that incorporates them 24 . These corrections 



involve bottom squarks and gluinos in the loops, but there are also possibly large corrections 



from stop and chargino loops. Both can be large since they grow as /itan (3 or A t /J,tan(3 24 



2a s nigfi tan (3 m t A t /itan/3 

37r max(m?,m? ,m? ) 87r 2 v 2 sin 2 (3 max(/i 2 , m^, m-^) 

Outside the decoupling limit, the reduced bb couplings of the h state are given in this case by 

9hbb ~ gAbb ~ tan /3(1 - A b ) (16) 

and can be thus significantly reduced or enhance^ depending on the sign of /x and, possibly, 
also At. This is exemplified in the left panel of Fig. |9j where the ratio Rbb = BR(/i — > 
bb) /BR(H SM — > bb) is shown as a function of the parameter \i tan (3 before the constraints of 
Table [T] The two branches in the histogram are due to the sbottom and stop contributions in 
which Rbb is increased or decreased depending on the sign of /i. 

A deviation of the partial h — )• 66 width will enter the total Higgs width, which is dominated 
by the bb channel, and change the R values for the different Higgs decay channels. A reduction 
of Rbb would thus lead to an enhancement of the 77 and the WW/ZZ branching fractions. 
Figure 10 shows the values of Rbb and _R 77 in which we observe a highly anti-correlated variation 



of the two ratios, with the exception of the cases where the opening of the decay h — > XX 
suppresses both branching fractions. The preliminary results from LHC and the Tevatron are 
overlayed. 

3.5.3 Sparticle contributions to the hgg and I177 vertices 

Scalar top quarks can alter significantly the gg — >■ h cross section as well as the h — > 77 decay 



width 25 . The current eigenstates ti,tR mix strongly, with a mixing angle oc mtX t , so that 



3 The other possible channels are strongly suppressed or ruled out. The co-annihilation with charginos, 
heavier neutralinos and staus is not effective as these particles need to be heavier than w 100 GeV and thus 
the mass difference with the LSP is too large. The annihilation through the A-pole needs Ma k, 2m x o < M/> 
and sizable tan /3 values, which is the anti-decoupling regime that is excluded as discussed above. Remains 
then the bulk region with staus exchanged in the t-channel in X1X1 ~* t + t~ (sbottoms are too heavy) which 
is difficult to enhance as the LSP is bino-like. 

4 These corrections also affect the Higgs production cross sections in the channels gg + bb — >■ <I>. However, in 
the cross sections times branching ratios for the t + t~ final states, they almost entirely cancel as they appear 
in both the production rate and the total Higgs decay width [57| . 
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Figure 9: (Left) Rbb values for a sample of pMSSM points as a function of the product of the //tan /3 
showing the reduction at large values of //tan/3. The reduction in a narrow strip at small values of \x 
is due to decays into XX- CRisht) the same as a function of At. 
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Figure 10: Rbb as a function of i? 77 , showing their anti-correlated variation; the points corresponding 
to a decrease of both ratios are due to an enhancement of invisible decays to light neutralinos. The 
values of i? 77 obtained by ATLAS+CMS and Rbb corresponding to the CMS and CDF+DO searches 
are overlayed as comparison. 



for large X t = A t — fi/ tan/3 values^] there is a lighter mass eigenstate t\ which can be much 
lighter than all other scalar quarks, <C M s . The coupling of the h boson to the ti states 
in the decoupling regime reads 



g hilk =cos2/3M| 



cos 2 6 t — -sL cos 26 1 



i 2 . 1 • oa v na ^rn t X t 

+ m. H — sin 2t) t m t X t , sin 26^ = — k 5- 

2 mi — mi 

tl *2 



(17) 



In the no-mixing scenario X t ~ 0, the coupling above is oc m\ and the scalar top contribution 
to the hgg amplitude is small, being damped by a factor 1/m? and interferes constructively 
with the top quark contribution to increase the gg — )■ h rate. However, since in the no-mixing 



5 On should assume X t values such that A t < 3Mg to avoid dangerous charge and colour breaking minima. 
In addition, if X t > \/6Ms, the radiative corrections to the h boson mass become small again and it would be 
difficult to attain the value M h « 126 GeV. 
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y/mti m t 2 ^ as ^° be very large for the h boson mass to reach a value Mh 



126 



scenario Ms 

GeV, the stop contribution to the hgg vertex, oc m|/M|, is very small. In the maximal mixing 
scenario, X t ~ \/6Ms, it is the last component of which dominates and becomes very 
large, oc — {mtXt/m^) 2 . However, in this case, the large contribution of a light stop to the 
hgg amplitude interferes destructively with the top quark contribution and the gg — >■ h cross 
section is suppressed. For ps 200 GeV and X t « 1 TeV, we obtain a factor of two smaller 
gg —7- h rate. In the case of sbottom squarks, the same situation may occur for large sbottom 
mixing Xf, = — /itan (5. However, for large value of Ms, it is more difficult to obtain a small 
enough m bi state to significantly affect the gg ^ h cross section. 

In the case of the /177 decay amplitude, there is the additional SM contribution of the W 
boson, which is in fact the dominant. Also, it has the opposite sign to that from the top quark 
and, hence, when stops are light and have a strong mixing, they will tend to increase the ^77 
amplitude. However, because the W contribution is by far the largest, the stop impact will be 
much more limited compared to the ggh case and we can expect to have only a ps 10% increase 
of the h — > 77 decay rate for m~ tl ps 200 GeV and X t ~ I TeV 25 -27 . Therefore, for light and 



strongly mixed stops, the cross section times branching ratio /x 77 is always smaller unity and 
relatively light stops do not entail an enhancement of the 77 yield. The sbottom contribution 
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Figure 11: i? 77 values for a sample of pMSSM points as a function of m^± (left) and (m r ±) (right). 
We impose Rbb > 0.9, to remove the effects due to the changes of the total width through the bb 
channel. 



to the /177 vertex is also very small, for the same reasons discussed above in the case of the hgg 
amplitude, and also because of its electric charge, — | compared to +| for stops. Other charged 
particles can also contribute to the h — >■ 77 rate 26 . The charged Higgs bosons have negligible 



contributions for m#± > 200 GeV. Charginos contribute to the hyy vertex and, because of 
their spin | nature, they contribution is only damped by powers of Mfe/m y ±. However, the 



/iXi 2X1^2 couplings are similar in nature to those of the LSP given in eq. (14) and cannot be 



strongly enhanced. As a result we expect contributions at most of the order of 10% even for 



mass values m x ± m 100 GeV (see Fig. 11). Charged sleptons have in general also little effect 
on the /i77 vertex, with the exception ofstaus |6]. These behave like the bottom squarks. At 
very large \i tan (3 values, the splitting between the two f states becomes significant and their 
couplings to the h boson large. Since f% can have a mass of the order of a few 100 GeV, without 
affecting the value of Mh, its contribution to the /177 amplitude may be significant for largte 



values of X T (see Fig. 11). 
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3.5.4 Impact of the LHC data 



Now, it is interesting to perform a first assessment of the compatibility of the LHC and Tevatron 
data with the MSSM and analyse the region of parameter favoured by the observed boson mass 
and rate pattern (see also [5j[62)). Despite the preliminary character of the results reported 
by the LHC collaborations and the limited statistical accuracy of these first results, the study 
is a template for future analyses. In this analysis, we computing the \ 2 probability on the 
observable of Table [ljfor each accepted pMSSM points. For the bb and r + r~ channels, in which 
no evidence has been obtained at the LHC, we add the channel contribution to the total x 2 only 
when their respective \x value exceeded 1.5 and the pMSSM point becomes increasingly less 
consistent to the limits reported by CMS. In order to investigate the sensitivity to the inputs, 
we also compare the results by including or not the bb, for which a tension exists between 
the CMS limit and Tevatron results, and the r + r~ rate. Figure 12 shows the region of the 
[X t , m 4 -J, [Xb, mgj and [Ma, tan f3] parameter space where pMSSM points are compatible with 
the input h boson mass and observed yields. In particular, we observe an almost complete 
suppression for low values of the sbottom mixing parameter X^. 
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Figure 12: Distributions of the pMSSM points in the [X t ,mjJ (left), rag ] (centre) and [Ma, 
tan j3] (right) parameter space. The black dots show the selected pMSSM points, those in light (dark) 
grey the same points compatible at 68% (90%) C.L. with the the Higgs constraints of Table [T] 

The distributions for some individual parameters which manifest a sensitivity are pre- 
sented in Figure 13, where each pMSSM point enters with a weight equal to its \ 2 probability. 
Points having a probability below 0.15 are not included. The probability weighted distri- 
butions obtained from this analysis are compared to the normalised frequency distribution 
for the same observables obtained for accepted points within the allowed mass region 122.5 
< Mjj < 127.5 GeV. We observe that some variables are significantly affected by the constraints 
applied. Not surprisingly, the observable which exhibits the largest effect is the product /i tan f3, 
for which the data favours large positive values, where the 77 branching fraction increases and 
the bb decreases as discussed above. On the contrary, it appears difficult to reconcile an en- 
hancement of both /i 77 and /x b g, as would be suggested by the central large value of [i bb - = 
1.97±0.72 recently reported by the Tevatron experiments [3]. Such an enhancement is not 
observed by the CMS collaboration and the issue is awaiting the first significant evidence of a 
boson signal in the bb final state at the LHC and the subsequent rate determination. The tan f3 
distribution is also shifted towards larger value as an effect of the Higgs mass and rate values. 
We also observe a significant suppression of pMSSM points with the pseudo-scalar A boson 
mass below ~450 GeV. This is due to the combined effect of the A — > t + t~ direct searches 
and B s — > fi + fi~ rate, which constrain the [Ma — tan f$\ plane to low tan value for light A 
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masses, by the shift to /xtan/3 from the Higgs rates disfavouring the low tan/3 region and by 
the suppression of the non- decoupling regime. 

In quantitative terms, we observe that 0.06 (0.50) of the selected pMSSM points are com- 
patible with the constraints given in Table [j] at the 68% (90%) confidence level. If we remove 
the constraint on the upper limit constraint on the bb and t + t~ rates, the fraction of points 
accepted at the 90% C.L. does not change significantly, at 0.56, but that at the 68% C.L. 
doubles to 0.12. On the contrary, if we replace the CMS upper limit for ji bb with the fi b b result 
of the Tevatron experiments for Mh = 125 GeV [3], the fraction of accepted points at 68% 
C.L. drops below 0.005. This highlights the tension which will be created in the pMSSM by a 
simultaneous excess in the 77 and bb channels, excess which cannot be adequately described 



in the pMSSM, as discussed above (see Figure 10) 
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Figure 13: The normalised distribution of the values of the /xtan/3 (upper left), tan/3 (upper right), 
Mi (lower left) and Ma (lower right) variables for the selected pMSSM points (dashed line) compared 
to the probability density function for the same variables obtained from the \ 2 probability using Mh, 
i? 77 and Rzz (continuous line). The normalised distributions reflect the biases induced by the Higgs 
constraints. 



4 Conclusions 

The implications of the new boson observation by the ATLAS and CMS collaborations for the 
phenomenological MSSM have been outlined. The study has been based on broad scans over the 
pMSSM parameter space where points have been preselected based on constraints from electro- 
weak and flavour physics, dark matter and searches at LEP2 and the LHC. Various scenarios 
for the stop mixing parameter X t (maximal, typical and zero-mixing) and representative values 
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of the soft SUSY-breaking scale M s (1, 2 and 3 TeV) have been confronted with the Higgs 
mass range compatible with LHC results, accounting for systematic uncertainties. In order to 
obtain Mh in the mass range 123 GeV< Mh < 129 GeV, large values of Ms and/or X t are 
required. In particular, the Mh constraints are sensitive to the value of the top quark mass 
for which the value extracted from the top quark pair production cross section has a more 
unambiguous definition but larger uncertainties. 

The various regimes of the pMSSM Higgs sector have been examined in the [M^tan/?] 
parameter. Of these regimes, only the decoupling regime, where the lighter h boson has 
almost SM-like properties and the heavier Higgs particles decouple from gauge bosons, and 
the SUSY regime survives all constraints. The anti-decoupling regime where the H state plays 
the role of the SM Higgs boson, the intense coupling regime in which there are three light 
states h, H and A, the vanishing coupling regime in which the h coupling to bottom quarks or 
gauge bosons are very strongly suppressed, and most of the intermediate coupling regime with 
relatively low Ma and tan /3 values, are excluded by the present data. In the SUSY regime light 
superparticles may affect the production and decay rates of the h boson. Light neutralinos 
may lead to invisible h boson decays, light stop and sbottom quarks affect the hbb couplings 
and the production cross section in the dominant gluon-gluon fusion mechanism, and light 
squarks, r-sleptons and charginos may affect the h — > 77 decay mode. 

We have confronted these possibilities with the recent LHC results and find that a significant 
fraction of pMSSM points in our scan compatible with them, including a possible enhancement 
of the 77 rate. Improved precision in the experimental measurements and sensitivity to the 
direct searches for the heavier Higgs bosons and supersymmetric particle partners at the LHC 
will provide the basis for clarifying the relation between the newly discovered scalar sector and 
physics beyond the Standard Model. 
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